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Plane Impact Response of PBX 9501 and Its Components below 2 GPa 


by 


Jerry J. Dick,∗ A. Richard Martinez, and Robert S. Hixson 


ABSTRACT 


The Hugoniot equation of state was measured for PBX 9501, the binder, and estane below 2 GPa using a 
light-gas gun facility. Time-resolved wave proﬁles were also obtained for these materials as well as for HMX 
single crystals in a state of uniaxial strain. The dynamic strength of PBX 9501 was measured at high strain rates 
in both compression and tension. In addition, elastic moduli were obtained from measurements of longitudinal 
and shear sound speeds. 


I. 
INTRODUCTION 


The problem of accidental detonation of explosive due 
to a mechanical insult is of continuing concern. The insult 
may not begin as a shock, but data from plane impact ex- 
periments provide equation of state data and the mechan- 
ical response at high strain rates in a well-deﬁned state of 
uniaxial strain. In these experiments we have measured 
impact response from 200 MPa to 2 GPa and impact ve- 
locities down to 120 m/s. There was essentially no existing 
data in this region for PBX 9501 or its components. The 
low end overlaps the range of stresses found in the spigot 
test, Steven test1, and other multidimensional tests.2 The 
results should enhance understanding of impact tests on 


PBX’s and serve as a basis for testing models used in com- 
putation of hazards problems. The experimental results 


are presented in considerable detail, so that workers in 


these areas will have a complete description of the experi- 
ments. 


Studying the components is of interest since there are 
models for computing the response of the composite ma- 


terial from the constitutive properties of the components. 


The binder is 5% of the PBX by weight and 7.2% by vol- 
ume. By weight the binder is half estane and half nitro- 


∗Electronic mail:jjd@lanl.gov 
∗Electronic mail:jjd@lanl.gov 


plasticizers. Actually, 0.1 wt % of an antioxidant is added 
as well. The response of estane itself is of interest to see if 
a good description can be obtained for it as an elastomer. 
There is interest in changes in molecular weight of estane 
over time. Migration of the nitroplasticizer is also of inter- 
est. 


II. 
EXPERIMENTAL TECHNIQUE 


Plane impact experiments were performed on PBX 9501 
and its components with shock strengths from 200 MPa to 
2 GPa (0.1 GPa=1 kbar=100 MPa=14.50 ksi) obtained 
using a light-gas gun facility. 
Argon was used for the 


propellant below about 400 m/s and helium above that. 


Time-resolved wave proﬁles of particle velocity were mea- 
sured using magnetic gauges in one of two conﬁgurations 


or using VISAR† laser interferometry. The experimental 


arrangements are depicted in Figs. 1 to 3. In the case of 
magnetic gauges, in-material wave proﬁles were measured; 


in the case of VISAR experiments the proﬁles were mea- 


sured at the interface between the material and a PMMA 


(polymethylmethacrylate)f window or vacuum. The laser 


beam is reﬂected from the mirror at the interface with spot 
size of about 0.5 mm at the center of the sample. The av- 


†VISAR 
stands 
for 
velocity 
interferometer 
system 
for 
any 
reﬂector. 
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erage target chamber temperature was 20.7 ± 2.3◦C. Den- 


sities were measured by water immersion. 


For low velocities the projectile body was made of brass 


with a Lexan polycarbonate nose and had a mass of about 


4.8 kg. For high velocities the projectile body was either 


Lexan or 2024 aluminum. The impactors for symmetric 


impact experiments with PBX 9501 were 50.8 mm in diam- 


eter. Three VISAR experiments 1053, 1054, and 1062 had 


impactors 50.8 mm in diameter ﬂush-mounted in an alu- 


minum projectile, but backed by PMMA 8 mm thick. The 


remaining experiments had impactors 57.2 mm in diameter 


and were backed by a microballoon foam or polyurethane 


foam with a density of about 0.48 g/cm3. 


The samples of PBX 9501 were 50.8 mm in diameter 
and the samples of estane and binder were 38 mm in di- 
ameter. The samples were bonded to the aluminum tar- 
get plate with epoxy, Sylgard 184 elastomer, or Aralhex 
polyurethane glue. For VISAR experiments the samples 
were placed in a 51 mm dia hole in the aluminum tar- 
get plate 6 mm thick. The impact face of the sample was 
nearly ﬂush with the front face of the plate; the impact 
face typically was 50 to 75 µm above the aluminum plate. 
For magnetic experiments the samples were mounted on 
the front of the plate so that the gauges were about 15 to 
25 mm above the plate. In this way the dynamic measure- 
ments were made before the shock wave reached the metal 
plate and disturbed the magnetic ﬁeld. 


A. 
PBX 9501 Samples 


With one exception, all PBX 9501 samples were from 
13.5 x 13.5 x 3.5-inch blocks hydrostatically pressed at S- 


Site from Holston lot 89C730-010 under the supervision 


of Manny Chavez. The drawback of this lot is that it was 
not possible to achieve densities within the desired range of 


1.828 to 1.838 g/cm3 by hydrostatic pressing. The average 


density of the samples was 1.827 ± 0.001 g/cm3. However, 
there was no immediate alternative available for hydrostat- 


ically pressed material. It was believed that for mechani- 


cal response this should not cause a signiﬁcant diﬀerence. 
Four possible variations are lot, pressing type, density, and 


age. The speciﬁcation density could be achieved by ram- 


die pressing. 
But ram-die pressing is believed to cause 


variation in density and particle size with radius and dis- 


tance from the ram. Furthermore, it is not the pressing 


technique of interest. 


B. 
Binder and Estane Samples 


The binder samples were fabricated at Building 360 at S- 


Site by Robert Montoya with supervision by John Kramer. 


A thin layer was poured onto a cookie sheet or into a muﬃn 


tin. Solvent was allowed to evaporate, and then another 


layer was poured. In this way samples up to 6 mm thick 


were made. The material is red, soft, and rubbery. Sam- 


ples 38 mm in diameter were cut out with a cork borer. The 
surfaces were ﬁnished by a surface grinder at room temper- 
ature. The samples were thickest in the center with about 
a 50 µm variation in thickness. Samples were bonded to 
a ﬂat PMMA window and then placed on a ﬂat surface. 
Then they were cast in epoxy so that a ﬂat impact face 
was achieved. The average and peak molecular weights of 
the estane in the binder were measured as 1.5 x 105 and 
1.2 x 105 g/mol by Wayne King. 


The estane samples were hot-pressed under vacuum to 
15 ksi at S-Site from lot 6812100SK13, DA95171-01 ﬂakes 
under the supervision of Cindy Sandoval. The ﬂakes are 
coated with talc to prevent sticking during storage. The 
pressing was 15 cm in diameter and 3.5 cm thick. Samples 
38 mm in diameter and up to 10 mm thick were machined 
at Hand Precision Machining by blowing liquid nitrogen 
over the surface to freeze it without freezing the main mass. 


Thickness varied by up to 110 µm on samples used in ex- 


periments. 
The estane sample disks are thickest in the 
center. Both average and peak molecular weights were .73 


x 105 g/mol in this form. 


For both binder and estane the PMMA window was 
bonded to the sample and the sample cast into a 51-mm 


dia hole in the target plate with epoxy. 


C. 
HMX Crystal Samples 


The HMX (cyclotetramethylene-tetranitramine) crystal 


slabs about 10 mm in lateral dimension were mounted on 
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Impactor 


PBX 
9501 


Inclined Gauges 
5 um Al in 25 um TFE sheets 


Inclined Gauges 
5 um Al in 25 um TFE sheets 
B, Magnetic Field 


EMV Stirrup Gauge 


Active Element 


Figure 1: Schema of the gun impact experiment 
with MMV instrumentation. Gauge positions are 
marked in red. The gauge leads and insulation 
are marked in blue. 


aluminum-coated, Kel-F∗ impact discs with a silicone elas- 


tomer. It was deemed important to not impact the crys- 
tals directly in order to avoid spurious eﬀects at the im- 
pact face. The crystal and PMMA window were bonded 
together with a silicone elastomer and surrounded with a 
casting of epoxy loaded with 40 vol % fused silica in order 
to reduce edge eﬀects. The epoxy loaded with silica is a 
good impedance match for explosives. 


The HMX crystal slabs were cut from larger crystals us- 
ing a low-speed, diamond saw. The crystals were polished 
using a series of polishing sheets consisting of alumina em- 
bedded in plastic. Soap and water were used as the carrier. 
The crystals were water-clear and inspected with a binoc- 
ular microscope at magniﬁcations up to 50x to make sure 
that the interiors were free of inclusions and other visible 
defects larger than about 2 µm. Crystals were grown by 
Howard Cady a number of years ago. 


D. 
Magnetic Gauging Experiments 


Magnetic gauges operate according to Maxwell’s third 


relation. The shock-induced motion of the gauge active 
element is perpendicular to the magnetic ﬁeld. Therefore, 


the induced voltage is ϵ = Bul, where B is the magnetic 


ﬁeld of about 0.750 kgauss, u is the particle velocity in- 
duced by the shock in mm/µs, and l is the length of the 


gauge active element in cm. The length is 0.6 to 1.0 cm. 


Two types of magnetic gauging experiments were per- 


∗Kel-F is polytriﬂuorochloroethylene manufactured by the 3M 
Company. 


Planar EMV Gauges 
4 um Al on 25 um TFE sheet 


Impactor 


Figure 2: Schema of the gun impact experiment 
with stirrup gauge instrumentation. Gauge po- 
sitions are marked in red. The gauge leads and 
insulation are marked in blue. 


formed on PBX 9501. Magnetic impulse and particle ve- 


locity (MIV) gauge packages consisting of 5 nested gauges 
of each type are placed side by side sandwiched between 
sheets of Teﬂon 25 µm thick (See Fig. 1.). 
The nested 
particle velocity gauges are 6 to 10 mm wide. 
The ac- 
tive elements were made of aluminum 5 µm thick. The 
gauge active elements are spaced 2 mm apart in the sand- 
wich so that when placed in the sample at an angle of 30◦ 


the gauges are successively 1 mm deeper into the explo- 
sive sample. Some experiments were performed with two 
sets of nested particle velocity (MMV, multiple magnetic 
velocity) gauges. 


Many of the earlier experiments had a nested set of tri- 
angular gauges designed to measure impulse in addition to 
the rectangular particle velocity gauges, MIV gauges (mag- 
netic impulse and velocity gauges). The intent was to use 
the two types of gauges in concert in order to perform La- 
grangian analysis on the data.3 In this way stress-strain 


curves could be obtained for loading and unloading at sev- 


eral depths in the sample. However, it has become clear 
that a robust, straightforward, Lagrangian analysis of that 


type requires a ﬁtting form for the data surfaces other than 


polynomials such as are used currently. 


The second type consists of two particle velocity gauges 


(EMV, electromagnetic velocity, stirrup) with 10-mm-wide 
active elements. These gauges are emplaced parallel to the 


shock impact face as in Fig. 2. One is placed on the impact 


face, and one is placed between two pieces of explosive of 
varying thicknesses. 
The gauge package consists of the 


aluminum gauge on one sheet of Teﬂon. 
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PZT Pin 


Laser Light 


Mirror 


PMMA Window 


Impactor 


Figure 3: Schema of the gun impact experiment 
with VISAR instrumentation. 


The samples of PBX 9501 were 50.8 mm in diameter. 


For magnetic experiments the samples were bonded to the 


face of the target plate and had free lateral surfaces. For 


stirrup gauge shots the PBX 9501 disks were bonded to a 


9 to 12 mm thick plastic disk in order to oﬀset the gauges 
from the aluminum target plate and avoid eddy currents in 
it due to shock-induced motion during the gauge recording 
time. 


E. 
VISAR Experiments 


For VISAR studies, PMMA (Mil. 
Spec. 
P-5425D, 
preshrunk) windows were used, Fig. 3. The mirrors were 
either vapor-plated aluminum or 13-µm-thick aluminum 
foil. 


The surface of the PMMA adjacent to the crystal was 
coated with a diﬀuse aluminum mirror for VISAR particle 
velocity measurements. 
Roughening the PMMA surface 
with 9-µm Al2O3 embedded in a plastic ﬁlm gave a satis- 
factory roughness for a diﬀuse reﬂector. Some repolishing 
can be done to get more light back into the receiving lens 


of the ﬁber optic probe. 
In some cases a 13 µm thick 


aluminum foil was used for the mirror instead of a vapor 


plated mirror. 


The samples of PBX 9501 were 50.8 mm in diameter 
and the samples of estane and binder were 38 mm in di- 


ameter. For VISAR experiments the samples were bonded 


to the PMMA window or aluminum foil mirror and glued 


in a 50 mm dia hole in the aluminum target plate. Two 


piezoelectric pins were emplaced on the impact face of the 
sample on opposite sides of a diameter. The time diﬀer- 


ence between shock arrival at these pins and the shock ar- 


rival signal recorded by the VISAR allowed measurement 


of transit time through the crystal and from that the wave 


velocities. It is important to account for the diﬀerence in 


signal travel times to the digitizing oscilloscopes for the 


two types of signals. 


The measurement system used was a push-pull, VISAR 


system.4 Some experiments were performed with two VIS- 


ARs and some with one VISAR. Using two VISARs with 


diﬀerent fringe constants removes ambiguity in determin- 


ing the particle-velocity jump at the shock when extra 


fringes must be added. This is important when initiation 


of the explosive sample occurs. The light was transported 


from the argon-ion laser to the target and thence to the in- 


terferometer table with ﬁber optics. The ﬁber optic probe 
on the target was obtained from Valyn International. It 
has a focal spot size of 0.6 mm. The light interference sig- 
nals were converted to voltage using photomultiplier tubes 
with rise and fall times of about 1 ns. The signals were 
recorded on digitizers at 1 or 2 ns sampling rates. Judging 
by the Lissajous ﬁgures generated by combining the sig- 
nals that are 90◦ out of phase, the overall response time of 
the system is about 3 ns at the shock-wave arrival time. 


F. 
Sound Speeds 


Sound wave velocities were measured with a pulse-echo 


apparatus, the same apparatus used by Stan Marsh and 


Bart Olinger at Ancho Canyon for many years. A burst 
of acoustic energy is emitted from a 5 MHz transducer. 


Round trip echo times are measured with a receiving trans- 


ducer. 
The moduli computed from the acoustic signals 
correspond to the limit of inﬁnitesimal strain and high 


strain rate. 
Materials measured were PBX 9501, PBX 


9501 binder, estane, HMX crystals, and Kel-F. Most of the 


materials studied were very dispersive. Therefore, signal 


averaging was used to measure signals from weak echoes. 
Sample thickness for sound speed measurement was typi- 


cally 1 to 1.5 mm for these materials. 
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III. 
SUMMARY OF EXPERIMENTAL 


RESULTS 


Wave proﬁles of particle velocity vs time were obtained 


for PBX 9501 at about 0.3, 0.6, 1.2, 1.6, and 2.0 GPa. For 


PBX 9501 there were no reliable Hugoniot data below 2 


GPa. Hugoniot equation of state data were obtained on 


PBX 9501, binder, estane, and HMX single crystals. For 


estane there were no data below 1 GPa. Furthermore, the 


intercept of the existing shock velocity vs particle velocity 


Hugoniot at zero particle velocity was not in agreement 


with the bulk sound velocity. No Hugoniot data existed 


for the binder, or for HMX single crystals below 35 GPa. 


In PBX 9501 multiple embedded magnetic gauge records 
were obtained at 0.3, 0.6 and 1.2 GPa with 5 particle ve- 
locity histories at 1.4 to 5.4 mm. These were symmetric 
impact experiments with PBX 9501 impactors. In this way 
the ﬁnal particle velocity is required to be half the accu- 
rately measured projectile velocity. This is a check on the 
accuracy of the MIV gauges. At 0.3 GPa 2 experiments 
with stirrup gauges at 5.4 mm depth were performed in or- 
der to check repeatability and agreement with the multiple 
embedded gauges that are in the sample on a plane slanted 
30◦ with respect to the shock plane. Similarly, one stirrup 
gauge experiment was performed at 0.6 GPa. At 0.3 and 
0.6 GPa stirrup gauge experiments were performed with 
8-mm thick PBX 9501 samples and VISAR experiments 
with 10 mm thick samples. In addition, wave proﬁles were 
obtained with VISAR instrumentation on samples 5 mm 
thick at 1.2, 1.6, 1.8 and 2.0 GPa. There was evidence of 


initiation at 1.6 to 2.0 GPa. 


For the binder and estane VISAR experiments were per- 


formed at about 200, 400, and 1000 MPa. From the ex- 


perimental results Hugoniot points and particle velocity 
histories were obtained. 


For HMX single crystals two VISAR experiments were 
performed to measure the strength of the elastic precursor 


in two crystal orientations. The two orientations were 110 


and 011 in the space group P21/n; the impact plane was 


parallel to those planes. Impact stress was about 1.4 GPa. 


Discussion of the several kinds of results on the diﬀerent 


materials follows. 


IV. 
HUGONIOT OF PBX 9501 


Hugoniot data for PBX 9501 are displayed in Fig. 4. A 


ﬁt to the Hugoniot data given in Table 1 and the data in 


Ref. 5 for densities 1.832 to 1.844 g/cm3 is 


US = 2.40 + 2.39up, 
0 < u < 0.9 mm/µs, 
(1) 


where US is shock velocity or wave velocity at half the 


maximum particle velocity, and up is ﬁnal particle velocity. 


The standard errors of the ﬁt parameters are 0.03, and 


0.07, respectively. The root mean squared (RMS) error of 
the shock velocity is 0.09 mm/µs. 
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Figure 4: PBX 9501 shock velocity vs particle 
Hugoniot data. The solid blue circle is the aver- 
age Hugoniot elastic limit point. 


The two data points below 600 MPa in Ref. 5 were not 


used as their validity is questionable. Those points were 
obtained using a streak camera and assumed a sharp shock. 


Our new wave proﬁle measurements show that this is not 


the case at such low longitudinal stresses. Rise times are 


on the order of a microsecond after several millimeters of 


wave propagation. 


The two spall experiments, 1082 and 1084, at 550 MPa 


were not used as the samples were only 3 mm thick. This 
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increases the uncertainty in wave velocity due to errors 


in measuring the transit time. 
Also, the wave is very 


unsteady in this region. 
The measured, apparent, wave 


speeds are low compared to other measurements at this 


stress level. If they are included the ﬁt is 


US = 2.37 + 2.44up, 
0 < u < 0.9 mm/µs. 
(2) 


The standard errors of the ﬁt parameters are 0.03, and 


0.08, respectively. The RMS error of the shock velocity is 


0.10 mm/µs. 


The wave proﬁles are unsteady up to 600 MPa for prop- 


agation distances up to about 8 mm. The Hugoniot jump 


conditions for conservation of mass and momentum ap- 
ply only to discontinuities or steady waves. 
The Hugo- 
niot curve is the locus of end states on the thermodynamic 
equilibrium surface for the material for shocks taking the 
material from the initial state to the ﬁnal state. Approxi- 
mate Hugoniot states were obtained by taking the transit 
time for the particle velocity on the wave front equal to 
half the ﬁnal particle velocity for the entire wave. It must 
be kept in mind that some of the Hugoniot data at or be- 
low 600 MPa are approximate because of the unsteadiness 
of the wave for the ﬁrst 8 mm of propagation distance; the 
wave proﬁle and rise time changes as it propagates in that 
distance. Using the particle velocity for the entire wave 
may systematically make the wave velocities near 0.3 GPa 
a bit high relative to the others. But it was decided to 
not overanalyze the data by trying to treat the unsteady, 
dispersive, wave front as two waves for Hugoniot points. 


The data are displayed as stress vs strain in Fig. 5. Vol- 


umetric strain is -(V − V0)/V0 = 1 − V/V0, where V0 is 
the speciﬁc volume in the initial state. The Hugoniot data 


represent the locus of end states for the compressive waves. 


The locus describes a curve on the thermodynamic surface 


of the material. Ten percent uniaxial strain is reached at 


about 1.8 GPa and 20% at about 7 GPa. 


Table 1 lists the measured Hugoniot data for PBX 9501. 


The measured shock velocity, US, is given in column seven. 


The Hugoniot state in PBX 9501 for the leading shock is 


computed using impedance matching techniques. For sym- 
metric impact experiments the particle velocity is half the 
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Figure 5: 
Longitudinal stress vs volumetric 
strain for PBX 9501. 
The volumetric strain is 
1-V/V0. 


impactor velocity. In other cases the particle velocity is 
computed from the intersection of the Rayleigh line for the 
sample and the reﬂected Hugoniot of the impactor origi- 
nating at the impactor velocity in the stress vs particle 
velocity plane. 


A 30◦ MMV experiment, 1075, by Rick Gustavsen on 
PBX 9501 with a density of 1.826 g/cm3 gave a particle 
velocity of 0.510 and a shock velocity of 3.64 for an impact 
stress of 3.39. 
This point agrees well with the existing 
wedge Hugoniot points in that region with higher PBX 
9501 densities. 


V. 
SOUND SPEEDS, ELASTIC MODULI, 


AND DENSITIES 


An isotropic material has just two independent elastic 


constants. Therefore, all elastic moduli can be calculated 
from two independent sound speeds, e.g., longitudinal CL 
and shear CS. The shear modulus G is ρ0C2 
S. The bulk 


modulus is 


K = ρ0C2 
0 = ρ0 


C2 
L − 4 


3C2 
S 


. 
(3) 
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Table 1: Results of Impact Experiments on PBX 9501 


Shot 
Experiment 
Impactor, 
Samplea 
Thick- 
num. 
type 
d 
ρ0 
ness 
uI 
US 
upb 
Pxb 
ϵxb 


(mm) 
(g/cm3) 
(mm) 
(mm/µs) 
(mm/µs) 
(mm/µs) 
(GPa) 


1115 
MMV 
9501, 3.02 
1.826 
1.33-5.27 
0.1325 
2.67 
0.0662 
0.323 
0.025 


1116 
MMV 
9501, 2.80 
1.825 
1.32-5.27 
0.1322 
2.51 
0.0661 
0.303 
0.026 


1071 
Stirrup 
Kel-F, 4.19 
1.827 
5.405 
0.133 
2.56 
0.066 
0.307 
0.026 


1072 
Stirrup 
Kel-F, 4.19 
1.828 
7.993 
0.134 
2.61 
0.065 
0.310 
0.025 


1061 
VISAR 
Kel-F, 5.89 
1.827 
10.018 
0.148 
2.59 
0.072 
0.341 
0.028 


1070 
Stirrup 
Kel-F, 4.19 
1.827 
5.415 
0.149 
2.57 
0.074 
0.345 
0.029 


1049 
MIV 
9501, 3.01 
1.827 
1.42-5.4 
0.233 
2.63 
0.116 
0.560 
0.044 
1117 
Stirrup 
z-qtz, 6.36 
1.826 
5.39 
0.1585 
2.58 
0.1236 
0.582 
0.048 
1073 
Stirrup 
Kel-F, 4.83 
1.827 
7.976 
0.248 
2.64 
0.120 
0.579 
0.045 
1060 
VISAR 
Kel-F, 5.89 
1.828 
10.018 
0.247 
2.73 
0.120 
0.599 
0.044 
1042 
MIV 
9501, 3.23 
1.835 
1.38-5.32 
0.233 
2.80 
0.116 
0.598 
0.042 
1082c 
VISAR 
z-qtz, 3.02 
1.827 
3.028 
0.1576 
2.44 
0.125 
0.556 
0.051 
1084c 
VISAR 
z-qtz, 3.03 
1.827 
3.01 
0.1530 
2.47 
0.121 
0.546 
0.049 
1058 
MIV 
9501, 2.99 
1.828 
1.37-5.33 
0.441 
2.93 
0.220 
1.18 
0.075 
1054 
VISAR 
z-qtz, 2.55 
1.826 
5.04 
0.300 
2.86 
0.230 
1.20 
0.080 
1053 
VISAR 
z-qtz, 2.57 
1.825 
5.016 
0.380 
2.97 
0.289 
1.56 
0.097 
1062 
VISAR 
z-qtz, 2.56 
1.828 
4.982 
0.381 
3.07 
0.287 
1.61 
0.094 
1124 
VISAR 
z-qtz, 3.010 
1.826 
10.015 
0.3754 
3.05 
0.2825 
1.57 
0.093 
1125 
VISAR 
z-qtz, 3.007 
1.828 
4.956 
0.4140 
3.20 
0.3080 
1.80 
0.096 
1081 
VISAR 
z-qtz, 6.37 
1.823 
4.538 
0.446 
3.22 
0.333 
1.96 
0.103 
1075d 
MMV 
— 
1.826 
— 
— 
3.64 
0.510 
3.39 
0.140 


aThe PBX 9501 samples were all hydrostatically pressed from Holston lot 89C730-010 except for the one used in shot 
1042. It was ram-pressed from another lot. 


bValues were calculated from the intersection of the Rayleigh line and the appropriate reﬂected impactor Hugoniot using 


the wave velocity at half maximum particle velocity, the measured impactor velocity and the known impactor Hugoniot. 


In cases of symmetric impact the particle velocity was taken as half the impactor velocity. In stirrup gauge shots the 


particle velocity was taken as that measured by the gauge on the impact face, and the wave velocity was taken from the 


transit time between gauges. For VISAR shots the state was computed from the measured wave velocity and ﬁnding 
the intersection of the Rayleigh line and the reﬂected impactor Hugoniot. 


cThe primary purpose of these shots was measurement of spall strength. The pulse from the impactor was 1 µs in 
duration, and the spall pullback occurred at the rear surface before the ﬁnal amplitude of the transmitted pulse was 


recorded at the rear free surface. Consequently, these factors along with the thin sample thickness and small transit 


time apparently reduce the accuracy of a Hugoniot point computed from this experiment. They appear to be outliers, 
and were not used in the Hugoniot ﬁt. 


dThese data were acquired by Rick Gustavsen using 30◦ MMV gauges. 
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Poisson’s ratio is 


ν = 0.5 



1 − 
1 


CL 
CS 


2 
− 1 



. 
(4) 


Young’s modulus is expressed as 


E = ρ0 
3C2 
L − 4C2 
S 


CL 
CS 


2 
− 1 
. 
(5) 


Longitudinal and shear wave velocities were measured 


for PBX 9501 binder, estane, HMX crystals and Kel-F. 


The results along with moduli are given in Table 2. Values 
for a type of copper are given for comparison. 


The moduli of polymers and PBX 9501 are rate depen- 
dent. From sound speeds the Young’s modulus for PBX 
9501 is 9.56 GPa. In contrast, a value of 0.77 GPa was 
measured in Instron tests.6 The ratio of high to low strain 
rate Young’s modulus is 12.4. The Young’s modulus for 
estane from Instron tensile data is about 7.6 MPa.7 In 
comparison the high strain rate value from sound speeds is 
340 MPa, 45 times greater. Noteworthy among the mod- 
uli are the low shear moduli and corresponding Poisson’s 
ratios near 0.5 for the binder and estane. These values are 
typical of rubber and elastomers. 


The longitudinal modulus, ρ0C2 
L, of PBX 9501 is 16.8 
GPa. For the two HMX crystal orientations the average 
value is 26.6 ± 0.3 GPa. In comparison, the c11 elastic 
constant of HMX is 20 GPa and the c33 is 21 GPa.8 


The sound speed values for PBX 9501 are the ones mea- 
sured by Bart Olinger some 20 years ago.9 His explosive 


was ram/die pressed. We tried to measure the longitudi- 


nal sound speed on some current material. PBX 9501 was 
very dispersive. A single echo was obtained with a sample 


1 mm thick. The shift in time of that echo was measured 


by adding another piece 1.5 mm thick. Three diﬀerent lots 


of PBX 9501 were tested in this manner. Discs 25 mm of 


the two thicknesses in diameter were machined from larger 
discs, typically 50 mm in diameter and 3 to 5 mm thick. 


Densities were measured on the larger discs at S-Site and 


on the small discs at TA-9. The results are shown in Ta- 


ble 3. The longitudinal sound speeds measured currently 


are all lower than Olinger’s. For comparison, the velocity 


of the leading foot of the elastic shock for 600 MPa im- 


pact stresses was about 2.92 for hydropressed Holston lot 


89C730-010. For one experiment from Holston lot 85G730- 


007, ram-pressed material with a density of 1.835 g/cm3 


the leading foot of the wave had a velocity of 3.00 mm/µs. 


The variation in the density measured at two diﬀerent 


times and places on the parent piece and on parts machined 


from it is larger than expected. Two 25 mm dia discs were 


machined from each of the three parent pieces. The varia- 


tion in density of the pairs of 25 mm dia pieces measured 


at TA-9 were 0.001, 0.004, and 0.003 g/cm3, respectively. 


VI. 
WAVE PROFILE MEASUREMENTS IN 
PBX 9501 


Wave proﬁles were measured in PBX 9501 at about 0.3, 
0.6, 1.2, 1.6, 1.8, and 2 GPa. Embedded, 30◦, nested mag- 
netic experiments were performed at 0.3, 0.6, and 1.2 GPa. 
The experiments at 0.3 GPa had a stirrup gauge emplaced 
on the center of the impact face. The measured proﬁles 
for the Holston lot 89C730-010 are shown in Figs. 6 to 9. 


At 0.31 GPa we see long rise times that evolve with 
propagation distance and a subtle two-wave structure de- 
veloping at the last two gauges (Figs. 6 and 
7). There 
is a change in slope just below 30 m/s. The ﬁrst wave is 


the elastic wave, and the change in slope deﬁnes the elastic 
limit. The inelastic wave is dispersive and unsteady. The 


unloading wave is recorded as well. The same behavior is 


noted at 0.6 GPa (Fig. 8). At 1.2 GPa the wave has sharp- 
ened up considerably, but there is an apparent shoulder at 


the end of the sharp rise (Fig. 9). At late times the agree- 


ment with half the projectile velocity is quite good at 1.2 


GPa, within 1% for the last 4 gauges. 
At 0.6 GPa the 


particle velocity at late times is 0.4 to 3.0% high. At 0.3 
GPa the ﬁnal state on the ﬁrst embedded gauge is 0.6% 


low in experiment 1116 and 4.7% low in experiment 1115. 


At subsequent gauges the ﬁnal state is not reached before 
the rarefaction arrives. 
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Table 2: Sound Speeds and Elastic Moduli 


Material 
Density 
CL 
CS 
C0 
ν 
K 
G 
E 


(g/cm3) 
(mm/µs) 
(mm/µs) 
(mm/µs) 
(GPa) 
(GPa) 
(GPa) 


PBX 9501 
1.82 
2.97 
1.39 
2.50 
0.36 
11.4 
3.52 
9.56 


binder 
1.27 
1.72 
0.28 
1.69 
0.49 
3.65 
0.100 
0.30 


estane 
1.19 
1.79 
0.31 
1.75 
0.48 
3.69 
0.115 
0.34 


HMX SX {110}n 
1.90 
3.76 


HMX SX {011}n 
1.90 
3.72 
Kel-F 
2.14 
1.89 
0.75 
1.68 
0.41 
6.04 
1.20 
3.39 
copper 
8.94 
4.75 
2.30 
3.94 
0.35 
138. 
47.3 
127. 


Table 3: PBX 9501 Sound Speeds 


PBX 9501 
pressing 
S-Sitea 
TA-9a 


type 
type 
density 
density 
CL 
(g/cm3) 
(g/cm3) 
(mm/µs) 


89C730-010 
hydro 
1.827 
1.821 
2.87 
unknown 
hydro 
1.834 
1.830 
2.86 


Blend 622-7 
ram 
— 
1.826 
2.78 


Olinger’s 
ram 
— 
1.82 
2.97 


a The S-Site densities are those of the parent pieces about 50 mm in diameter and about 5 mm thick, whereas the TA-9 


densities are of the 25 mm diameter by 1.5 mm thick pieces used in the sound speed measurements that were machined 
from the parent pieces. 
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Figure 6: In-material MMV gauge records in 
PBX 9501 at 303 MPa in a symmetric impact 
experiment, 1116. 
Gauge positions are 1.32 to 
5.27 mm in from impact face. Half the impact 
velocity is 66.2 m/s. The middle part of the sec- 
ond record at 2.32 mm was lost. The impact-face 
gauge leads were potted in epoxy ﬂush with the 
impact face. The hydrodynamic ﬂow near the lat- 
eral edge apparently caused gauge-lead spreading 
and anomalous increases in voltage level. There- 
fore, only the ﬁrst part of the impact face record 
is displayed. 
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Figure 7: In-material MMV gauge records in 
PBX 9501 at 323 MPa in a symmetric impact 
experiment, 1115. 
Gauge positions are 1.33 to 
5.27 mm in from impact face. Half the impact 
velocity is 66.1 m/s. 
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Figure 8: 
In-material MIV gauge records in 
PBX 9501 at 560 MPa in a symmetric impact 
experiment, 1049. 
Gauge positions are 1.38 to 
5.32 mm in from impact face. Half the impact 
velocity is 116 m/s. 
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Figure 9: 
In-material MIV gauge records in 
PBX 9501 at 1.18 GPa in a symmetric impact 
experiment, 1058. 
Gauge positions are 1.37 to 
5.33 mm in from impact face. Half the impact 
velocity is 220 m/s. The origin of time is calcu- 
lated as the distance to the ﬁrst gauge divided by 
the longitudinal sound speed. 
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A. 
Wave Proﬁle Measurements in PBX 9501 


near 310 MPa 


In addition to the MMV gauging measurements, three 


stirrup gauge measurements were made (two at 5.4 mm 


and one at 8 mm) as well as a VISAR measurement with 


a 10 mm thick sample. 
The two stirrup gauges at 5.4 


mm and at 307 and 345 MPa agree very well with the 


MMV record at 5.27 mm from experiment 1116 at 303 


MPa (Fig. 10). This shows that the proﬁle measurements 


are repeatable on PBX 9501 at the same density, from the 


same lot, and pressed in the same way. This also demon- 


strates that gauges emplaced at 30◦ to the impact plane 


agree well with gauges emplaced parallel to the impact 
plane and the wave front. 


The VISAR experiment measures the particle velocity 
history at the interface between the explosive and the 
PMMA mirror. In addition, the experiments did not have 
exactly the same impact stress. 
In order to display all 
proﬁles on the same plot, it was decided to normalize the 
measured particle velocities by the expected ﬁnal particle 
velocities. For symmetric impact the ﬁnal particle velocity 
is half the projectile velocity; for VISAR and experiments 
with two stirrup gauges it was the value calculated from 
the known Hugoniots including the Hugoniot of PBX 9501 
using Joe Fritz’s MACRAME characteristics code. 


A summary of experiments near 310 MPa for the same 
lot of PBX 9501 is displayed in Fig. 11. A front-surface 
impact record is shown from experiment 1072. 
The in- 
terferometric record for the VISAR experiment at 10 mm 


was only half a fringe for the fringe constant used. This 


short of a record resulted in poor accuracy for the particle 
velocity proﬁle obtained. The ﬁnal level was about 10% 


low. Therefore the record was scaled by 1.10 in order to 


bring it in line with the other records and the ﬁnal particle 
velocity calculated by MACRAME. 


B. 
Wave Proﬁle Measurements in PBX 9501 
near 580 MPa 


The results of experiment 1049 at 560 MPa are displayed 
in Fig. 8. The rise times are sharper than at 310 MPa but 
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Figure 10: In-material magnetic gauge records 
in PBX 9501 near 310 MPa. It shows the relative 
particle velocity proﬁles for two stirrup gauge ex- 
periments, 1070 and 1071, at 5.4 mm depth com- 
pared to the MMV gauge proﬁle at 5.27 mm from 
experiment 1116 (magenta). 
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Figure 11: Experimental wave proﬁles in PBX 
9501 near 310 MPa for lot 89C730-010 hydrostat- 
ically pressed to a density of 1.827 ± 0.001 g/cm3. 
The two orange proﬁles are from stirrup gauges 
for experiment 1072. The last record is from a 
VISAR experiment 1061. 
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the wave is still unsteady, the proﬁle changing as it prop- 


agates. A stirrup gauge experiment 1117 was performed 


at about 5.37 mm for comparison with the last gauge in 


experiment 1049 that had the gauges emplaced at 30◦ to 


the impact plane at 5.39 mm. The comparison is shown 


in Fig. 12. The agreement in wave front proﬁles is good, 


similar to the results at 310 MPa. Both types of gauge 


emplacements perform well. 
The slanted gauges appear 


to show a sharper change in slope between the elastic and 


inelastic waves at a relative particle velocity of about 0.25. 


A summary of all records near 580 MPa plotted as relative 


particle velocity is shown in Fig. 13. 
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Figure 12: Experimental wave proﬁle records in 
PBX 9501 at 560 to 580 MPa. 
Comparison of 
wave proﬁles at 5.4 mm depth from a magnetic 
gauge embedded at 30◦ to the impact face (ma- 
genta curve) and one embedded parallel to the 
impact face (black curve). 


C. 
Wave Proﬁle Measurements in PBX 9501 at 
1.2 GPa 


The comparison of the nested, embedded, magnetic 
gauge records from experiment 1058 at depths from 1.37 


to 5.33 mm are compared to the VISAR record for exper- 


iment 1054 at 1.20 GPa with a sample 5.04 mm thick in 


Fig. 14. The rise times agree quite well, and the VISAR 


record conﬁrms the hesitation in the wave front rise at 92 
to 95% of the ﬁnal state. 


The z-cut quartz impactor in the VISAR experiment was 
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Figure 13: Experimental wave proﬁles in PBX 
9501 at 560 to 600 MPa. The two orange proﬁles 
are for experiment 1073. The last record is from 
a VISAR experiment 1060. 


only 2.55 mm thick. The shock velocity in the impactor 
was 6.46 mm/µs; therefore the rarefaction overtook the 
wave front before the ﬁnal state was reached after 5 mm 
of wave propagation in PBX 9501. 


The question arises whether the elastic wave is over- 
driven by 1.2 GPa. The velocity measured for the begin- 
ning of the wave rise is 2.94 mm/µs in experiment 1054 
and 2.99 mm/µs in experiment 1058. Thus the longitudi- 
nal elastic wave is on the verge of being overdriven at 1.2 
GPa. Results of experiments by Bob Graham and G. T. 
Holman at Sandia using polyvinylidene ﬂuoride (PVDF) 
gauges show evidence of a two-wave structure in the deriva- 


tive of the particle velocity at 1.2 GPa in PBX 9501. For 


comparison, the measured speed of the leading disturbance 
is 3.09 ± 0.06 mm/µs at 1.58 GPa, 3.29 mm/µs at 1.80 


GPa, and 3.27 mm/µs at 1.96 GPa. 


D. 
Wave Proﬁle Comparison in PBX 9501 at 0.3 
to 1.2 GPa 


A comparison of the wave front rise times and proﬁles 
from 0.3 to 1.2 GPa after about 5.3 mm of wave propa- 


gation is displayed in Fig. 15. These are records from the 


multiple nested magnetic gauge experiments 1116, 1049, 
and 1054. It shows the dramatic sharpening of the wave 
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Figure 14: Experimental wave proﬁles in PBX 
9501 at 1.2 GPa. The nested magnetic gauges are 
at 1.37 to 5.33 mm depth. The VISAR record is 
at 5.04 mm. 


front as the impact stress increases. It is typical of many 
materials that the rise time decreases and the strain rate 
increases as the impact stress increases. 


It should be mentioned that the time of impact (time 
zero) was measured with a stirrup gauge at 0.3 GPa. For 
the other two experiments it was computed from the posi- 
tion of the ﬁrst gauge divided by the appropriate longitu- 
dinal sound speed for deﬁning the time of ﬁrst motion at 
the ﬁrst gauge. 


E. 
Elastic Wave Strength in PBX 9501 


The elastic wave is discernible in experiments near 300 
and 600 MPa. In some measurements the end of the elastic 


wave is delineated by a subtle two-wave structure, in oth- 


ers by a change in slope of the wave proﬁle. The measured 
Hugoniot elastic limits (yield point in uniaxial compressive 


strain at high strain rates) are listed in Table 4. Corre- 


sponding yield stresses in uniaxial stress Y are calculated 


from10 


Y = 
1 − 2ν 


1 − ν 


σHEL 
x 
. 
(6) 


For samples from lot 89C730-010 the average Hugoniot 
elastic limit is 123 ± 8 MPa for impact stresses of about 


310 MPa and 150 ± 15 MPa for impact stresses of about 
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Figure 15: Experimental wave proﬁles in PBX 
9501 at 0.30, 0.56, and 1.20 GPa. The magnetic 
gauges are at about 5.3 mm depth. The proﬁles 
are from experiments 1116, 1049, and 1058. 


580 MPa. The corresponding yield stresses are 54 and 66 
MPa. The results for experiment 1042 at 600 MPa using 
material that was ram-pressed from another lot has the 
largest Hugoniot elastic limit by a signiﬁcant amount. It 
had a density of 1.835 g/cm3. 


For comparison, a Hugoniot elastic limit of 140 MPa 
was measured in Composition B.11, 12 A very diﬀerent, 
ramp-like precursor was observed for PBX 9502.13 The 
longitudinal stress value at the end of the ramp was 73 


MPa. Chhabildas and Kipp measured the maximum shear 


stress in PBX 9404 using pressure-shear impact loading.14 


Their result was 40 MPa. The yield stress is twice that, 


80 MPa. From Vickers hardness tests a value of 66 MPa 


was given for the ﬂow stress for HMX by Mohan, Bhasu, 
and Field15 giving a value of 132 MPa for Y. 


It is interesting to note that the threshold velocity for 
explosion in the spigot test or the Steven test corresponds 


to about 100 MPa in PBX 9501 in the center after im- 


pact, close to the Hugoniot elastic limit. The threshold 
velocity may be correlated to the threshold for inelastic, 


irreversible, dissipative deformation. 
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Table 4: PBX 9501 Hugoniot Elastic Limit and Yield Stress 


Hugoniot elastic limit 


Shot 
Particle 
Yield 
Total rise 
num.a 
velocity 
Stress 
stress 
U el 
0 
U el 
1/2 
time 


(m/s) 
(MPa) 
(MPa) 
(mm/µs) 
(mm/µs) 
(ns) 


1071b 
26 
126 
55 
2.73 
2.66 
160 


1072 
22 
113 
50 
2.90 
2.82 
160 


1070 
25 
129 
57 
2.87 
2.82 
150 


1049 
27 
140 
62 
2.90 
2.85 
90 


1073 
31 
161 
71 
2.94 
2.84 
140 


1042c 
41 
219 
96 
3.00 
2.91 
75 


aThe ﬁrst three experiments had impact stresses around 300 MPa, and the last three around 600 MPa. Experiments 
1072 and 1073 had samples 8 mm thick. 


bIn this experiment the rise time on the stirrup gauge on the impact face had a rise time of 120 ns. 


cThe sample in this experiment was ram-pressed from another lot. 


F. 
PBX 9501 Strength vs Strain Rate 


Using Instron and Hopkinson bar instrumentation Los 
Alamos National Laboratory workers have measured the 
ultimate compressive strength on the same lot of PBX 9501 
pressed to the same density by the same technique.6, 16 


The ultimate compressive strength is the maximum com- 
pressive stress achieved. 
Yield stresses are proportional 
to and less than the ultimate compressive strengths. The 
results are displayed in Fig. 16. 


Earlier measurements of PBX 9501 ultimate compres- 
sive strength are listed in Ref. 5 for several temperatures. 


There is strong temperature dependence of strength. The 


strain rate is not speciﬁed but may have been 0.00083 s−1. 


The density is 1.844 g/cm3. At 24◦C the ultimate compres- 


sive strength is given as 7.86 MPa. This value is slightly 
lower than values for the current material. In other earlier 


work on another lot of PBX 9501 reported in Dobratz’s 


compendium,17 the ultimate compressive strength in In- 
stron tests is about 6 MPa at a strain rate of 10−5 s−1, 33 


MPa at 750 s−1, and 43 MPa at 103 s−1. 


There is signiﬁcant strain-rate dependence of the yield 


stress. The state of strain may also be important. Mechan- 


ical failure in PBX 9501 and other PBXs at low strain rates 
has been studied in detail by Wiegand.18 


G. 
Wave Proﬁles in PBX 9501 from 1.6 to 2 GPa 


There is evidence of initiation for impact stresses from 
1.6 to 2 GPa. 
The experiments were 1053 and 1062 at 
1.56 and 1.61 GPa and 5 mm sample thickness, 1124 at 
1.57 GPa and 10.0 mm, 1125 at 1.80 GPa and 4.96 mm, 
and 1081 at 1.96 GPa and 4.54 mm. A summary of the 
proﬁles obtained is displayed in Fig. 17. 
All impactors 


were z-cut quartz. They were 57.2 mm in diameter and 
partially supported laterally with epoxy on all shots ex- 


cept 1053 and 1062 where they were 50.8 mm in diameter 


and ﬂush-mounted in the aluminum projectile with lateral 
aluminum support. The experiments were one-dimensional 


(uniaxial strain) at the mirror center until in the neighbor- 


hood of 5 µs. The quartz impactors were backed by foam 


with a density of about 0.48 g/cm3 and a thickness of 5 
or 6 mm, except for experiments 1053 and 1062 that were 


backed by 8 mm of PMMA. The results represent delayed 


initiation. As such, they should provide good test data for 
numerical modelling of hazards tests that are performed 


at lower stresses with longer delay times for initiation and 
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Figure 
16: 
Ultimate compressive strengths 
and yield stresses vs strain rate for PBX 9501. 
The square blue points are ultimate compres- 
sive strengths from Instron tests at about 27◦C. 
Error bars are less than 0.5 MPa, smaller than 
the points. Specimens were 9.5 mm diameter by 
19 mm long. 
The ends were lubricated. 
The 
red square is an ultimate compressive strength 
from Hopkinson bar tests at about 16◦C. Sam- 
ples were 6.35 mm in diameter and 3.18 mm 
long. 
The black circles are yield stresses com- 
puted from Hugoniot elastic limits in plate im- 
pact tests at about 23◦C. All samples are from 
the same lot and hydropressed to the same den- 
sity 1.827 g/cm3. 


for modelling of the hazard scenarios themselves. 


At the highest impact stress of 1.96 GPa, experiment 
1081, (upper black curve in Fig. 17) strong initiation oc- 
curs during a shock pulse of about 2 µs duration; the 
impactor was 6.37 mm thick. The initial particle veloc- 


ity at the PMMA window is calculated to be 402 m/s by 


MACRAME with an interface longitudinal stress of 1.54 


GPa. The exothermic decomposition leading to initiation 


is made evident by the rising particle velocity and pressure 
in the wave proﬁle from the VISAR record for a sample 


4.538 mm thick. Strong decomposition begins about 0.8 


µs after shock arrival, well before the rarefaction arrives 
from the back of the impactor. The Pop plots in LASL 


Explosive Property Data5 give a run distance of 21.2 mm 


and time to detonation of 6.45 µs for a density of 1.833 


g/cm3. For a density of 1.844 g/cm3 it gives a run dis- 


tance of 36.9 mm and a time to detonation of 10.5 µs. The 
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Figure 17: Experimental wave proﬁle records 
in PBX 9501 at 1.6 to 2 GPa at the PBX 
9501/PMMA interface using VISAR instrumen- 
tation. The PBX 9501 samples are about 4.5 to 
5 mm thick except for one at 10 mm. 


run distance for the density 1.823 g/cm3 of the sample used 
in experiment 1081 will be shorter than, but close to that 
for 1.833 g/cm3 material. 


Using a theoretical maximum density 1.860 g/cm3 the 
porosity is 2.0% for 1.823 g/cm3, 1.8% for 1.827 g/cm3, 
1.5% for 1.833 g/cm3, and 0.9% for 1.844 g/cm3. 
The 
1.844 g/cm3 material is close to dead pressed; that is, the 
concentration of hot spots is so low that it may be initiat- 
ing more like a homogeneous explosive. The run distance 
and time to detonation for 1.827 g/cm3 material should be 
somewhat shorter than for 1.833 g/cm3 PBX 9501. 


The next lower experiment, 1125, had an impact stress of 


1.80 GPa (red curve in Fig. 17). Its impactor was 3.00 mm 


thick. The initial particle velocity at the PMMA window is 
calculated to be 363 m/s by MACRAME with an interface 


longitudinal stress of 1.42 GPa. There is little evidence of 


initiation until late times with this lower impact stress and 
shorter pulse. The reduction in stress by the rarefaction 


has a signiﬁcant eﬀect on the initiation progress. 


Two nearly identical experiments, 1053 and 1062, were 


performed at about 1.60 GPa (blue and orange curves in 


Fig. 17). 
The impactor velocities were 0.380 and 0.381 
mm/µs, respectively. The impactors were 2.56 mm thick 
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backed by 8 mm of PMMA. The initial particle velocity 


at the PMMA window is calculated to be 333 m/s by 


MACRAME with an interface longitudinal stress of 1.30 


GPa. The two interface proﬁles agree well with each other 


to 4 µs and within 10% to 6.6 µs. The proﬁles show ev- 


idence of late initiation as the particle velocity increases. 


Another experiment 1124 was executed at 1.57 GPa (lower 


black curve in Fig. 17). The sample was 10 mm thick, and 


the impactor was 3 mm thick backed by foam. The pur- 


pose was to see the eﬀect of a longer distance of propa- 


gation on the initiation. The proﬁle shows that the rar- 


efaction has overtaken the shock front and attenuated the 


wave. Furthermore, the levels following the leading part 


of the rarefaction are lower than on the experiments with 
samples 5 mm thick. 
This is due to the lower particle 
velocities and stresses generated by the rarefaction from 
the interface between the quartz and the backing material. 
The PMMA backing material in the 5 mm experiments 
has higher impedance than the foam, and leads to higher 
particle velocities and stresses behind the leading rarefac- 
tion. The fact that the particle velocity begins to increase 
at later times with the 10 mm thick sample than the 5 mm 
ones could mean that the initiation begins at the impact 
face, and has farther to travel to reach PBX 9501/PMMA 
measurement interface. Also, the lower stresses and par- 
ticle velocities behind the lead rarefaction could cause a 
longer induction time. 


In most of these experimental proﬁles the proﬁles are 
of interest to late times. In hazard situations the onset 


of initiation of explosion or detonation may occur after a 
considerable delay from the time of the initial mechanical 


stimulus. The experiments presented here are completely 


one-dimensional at the measurement interface center for 
about 5 µs after impact. Two-dimensional numerical cal- 


culations would be necessary to determine the deviation 


from one dimensionality out to 9 µs. Furthermore, the ef- 


fect of two-dimensional ﬂow on the VISAR signal is not 


known. However, we had a u-t proﬁle at 1.2 GPa where 
we were fairly certain that initiation was not taking place. 


If our memory is correct there was residual explosive in 


the target chamber after the shot, no soot. Based on these 
considerations it seemed worthwhile to compare the exper- 


iments to one-dimensional calculations of the u-t proﬁles 


at the measurement interface using the MACRAME char- 


acteristics code. 


The comparison for experiment 1054 at 1.2 GPa is dis- 


played in Fig. 18. There is rough agreement out to 8 µs. 


The experimental proﬁle does not reach the maximum par- 


ticle velocity because at 1.2 GPa and 5 mm propagation 


distance, PBX 9501 does not reach the ﬁnal state until 


several hundred nanoseconds behind the main compressive 


wave (see Fig. 14). The rarefaction overtakes it before the 


ﬁnal state is reached. MACRAME is a hydrodynamic char- 


acteristics code. This means the shocks are discontinuities 


and the lead rarefaction moves at bulk sound speed. In a 
composite material with strength such as PBX 9501 the 
shock is dispersed, and the lead rarefaction moves at lon- 
gitudinal sound speed, faster than the bulk sound speed. 
This explains some of the diﬀerences between the com- 
puted and measured proﬁles. Still, the agreement in par- 
ticle velocity levels behind the leading rarefaction is good 
out to almost 8 µs. The jump in particle velocity in both 
proﬁles at about 7.5 µs is noteworthy. This jump is caused 
by the shock from the aluminum projectile body that sup- 
ports the 8mm thick PMMA disc that is behind the quartz 
impactor. This shock will appear for the impactors backed 
by PMMA but not for those backed by foam because of 
higher PMMA wave speeds. The overall experience from 
the VISAR experiments is that the records may have va- 
lidity up to about 9 µs, at least for z-cut quartz impactors 
no more than 3 mm thick. 


All the materials except z-cut quartz have wave speeds 


around 3 mm/µs. So, for the samples with 25.4 mm ra- 
dius the lateral rarefactions should not reach the center 


until after 8 µs. 
The z-cut quartz has velocities over 6 


mm/µs. So, if it were isolated the rarefactions could reach 


the center of it in about 4 µs. However, for thin impactors 


in the experiments there are reverberations in it from the 
surrounding PBX 9501 and foam. These cause drag on the 


lateral eﬀects and minimize them until a much later time 


after impact. 


Based on these observations at 1.2 GPa it is worthwhile 


to make the same sort of comparison for the experiments 
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Figure 18: Comparison of the wave proﬁle for 
experiment 1054 in PBX 9501 at 1.2 GPa at the 
PBX 9501/PMMA interface with a calculation by 
the MACRAME characteristics code. The PBX 
9501 sample is 5.04 mm thick. The z-cut quartz 
impactor was 2.55 mm thick and backed by 8 mm 
of PMMA. 


at higher impact stresses. The result for experiment 1053 
at 1.56 GPa and 5.02 mm is shown in Fig. 19. The ini- 
tiation becomes discernible at about 3 µs and accelerates 
at about 5 µs. This is within the one-dimensional time 
for the experiment. Without initiation the calculational 
results imply that the interface is held at 0.6 to 0.8 GPa 
when the quartz impactor is backed by PMMA. 


For experiment 1124 at 1.57 GPa and 10.02 mm the 
comparison is shown in Fig. 20. This comparison shows 


clearly the attenuation of the thin shock pulse in the ex- 


periment with wave propagation distance. The wave front 
is attenuated by about 12%. The impactor is backed by 


low impedance foam. 
The calculated proﬁle shows that 


the particle velocity drops continually with foam backing 


as waves reverberate across the impactor, in contrast to 


PMMA backing. This same behavior of more rapid and 
continuous reduction in pressure will have taken place at 


the PBX 9501 impact face. At about 4 µs the calculated 


interface pressure drops to about 0.65 GPa and at about 
5 µs to about 0.4 GPa. The experimental proﬁle begins to 


stay at higher particle velocities than the calculated pro- 


ﬁle at about 5 µs and shows strong apparent initiation at 
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Figure 19: 
Comparison of the experimental 
wave proﬁle 1053 in PBX 9501 at 1.56 GPa at the 
PBX 9501/PMMA interface with a calculation by 
the MACRAME characteristics code. The PBX 
9501 sample is 5.02 mm thick. The z-cut quartz 
impactor is 2.57 mm thick and backed by 8 mm 
of PMMA. 


about 6 µs. 


The comparison for experiment 1125 at 1.80 GPa is dis- 
played in Fig. 21. The particle velocity and pressure drop 
rapidly as waves in the impactor reverberate between the 
PBX 9501 and the foam. There may be evidence of ini- 
tiation beginning at 2.5 µs and clear evidence after 3 µs. 
The computed interface pressure drops to about 0.72 GPa 
at 2.5 µs and to about 0.45 GPa at 3.5 µs. At 5.2 µs it 
drops to 0.2 GPa; this is from the reverberation of the ini- 
tial shock in the sample. At 6.2 µs the pressure increases 


to 0.35 GPa; this is from the reverberation of the initial 
rarefaction in the sample. The bump and rapid rise in par- 


ticle velocity in the measured proﬁle begins at about this 


time. 


In order to make a more convincing case for the pro- 


ﬁles, it was decided to perform an experiment with an 
inert sample holding the experimental geometry constant. 


Kel-F (polychlorotriﬂuoroethylene) plastic was used for the 


sample. It has a Px − up Hugoniot very close to that of 


PBX 9501. The impactor was z-cut quartz 3.015 mm thick 


backed by 5 mm of foam. 
The Kel-F sample was 5.01 
mm thick. These are a close match to experiment 1125 
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Figure 20: 
Comparison of the experimental 
wave proﬁle 1124 in PBX 9501 at 1.57 GPa at the 
PBX 9501/PMMA interface with a calculation by 
the MACRAME characteristics code. The PBX 
9501 sample is 10.02 mm thick. The z-cut quartz 
impactor was 3.010 mm thick and backed by 5 
mm of foam. 
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Figure 21: 
Comparison of the experimental 
wave proﬁle 1125 in PBX 9501 at 1.80 GPa at the 
PBX 9501/PMMA interface with a calculation by 
the MACRAME characteristics code. The PBX 
9501 sample is 4.96 mm thick. The z-cut quartz 
impactor was 3.007 mm thick and backed by 5 
mm of foam. 


with PBX 9501. The impactor velocity was 0.403 mm/µs 


for an impact stress of 1.66 GPa and particle velocity of 


0.307 mm/µs. The experimental results and comparison 


to MACRAME calculation are presented in Fig. 22. The 


experimental proﬁle for the Kel-F polymer is dispersed at 


the top of the compressive wave. This allows the rarefac- 


tion to overtake and attenuate it. The agreement for the 


thin pulse is good out to 7.5 µs. The MACRAME x-t re- 


sults show the arrival of the reverberation of the rarefaction 


wave in the sample by this time. Wave interactions at late 


times are complex and may get smoothed out. This sort 


of behavior is manifest in Fig. 18 at 5.5 µs. Wave speeds 


in Kel-F are somewhat slower than in PBX 9501 so that 


reverberation eﬀects appear at later times in Kel-F. Over- 
all one can conclude that the longitudinal particle velocity 
at the rear of the sample is not signiﬁcantly aﬀected by 
lateral ﬂow in the ﬁrst 8 µs after impact when thin z-cut 
quartz impactors are used. 
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Figure 22: 
Comparison of the experimental 
wave proﬁle 1136 in Kel-F at 1.66 GPa at the 
Kel-F/PMMA interface with a calculation by 
the MACRAME characteristics code. The Kel-F 
sample is 5.01 mm thick. The z-cut quartz im- 
pactor was 3.015 mm thick and backed by 5 mm 
of foam. 


Fig. 23 shows the comparison between the MACRAME 


calculation and the data at 1.96 GPa for experiment 1081. 


This experiment had a z-cut quartz impactor 6.37 mm 
thick providing a shock pulse over twice as long as for the 


other experiments in this set. At this stress the initiation is 
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Figure 23: Comparison of the wave proﬁle for 
experiment 1081 in PBX 9501 at 1.96 GPa at the 
PBX 9501/PMMA interface with a calculation by 
the MACRAME characteristics code. The PBX 
9501 sample is 4.54 mm thick. Measured sample 
density was 1.823 g/cm3. The z-cut quartz im- 
pactor was 6.37 mm thick and backed by 6 mm 
of foam. 


prompt within 0.8 µs of the shock. The successive particle 
velocity decreases in the calculation are seen in the mea- 
sured proﬁle superimposed on the particle velocity increase 
with time due to initiation. In the calculation the initial 
shock pressure at the interface is 1.54 GPa. After 3 µs it 
drops to 0.78 GPa, after 4 µs to 0.64 GPa, and after 5.5 µs 
to 0.35 GPa. The behavior of the measured proﬁle after 6 
µs is peculiar in dropping down and then going up again. 
The reverberation of the rarefaction through the sample 


arrives at the interface at about 6.3 µs in the calculation. 
While there is a small dip in the MACRAME proﬁle at 


about this time, the graph is terminated at 6 µs due to the 


increased eﬀect of lateral rarefactions in the thick quartz 
impactor. 


H. 
PBX 9501 Rise Times and Strain Rates 


The rise times of the proﬁles in the experiments in PBX 


9501 were measured. The rise time was taken as the rise 
time from 10% to 90% of the total signal rise. Since the 


approach to ﬁnal state is so gradual at low stresses, it is 


diﬃcult to pick oﬀ the precise time at times closer to the 


ﬁnal state than 90%. The values are recorded in Table 5 


and displayed in Fig. 24. One notes that at 0.3 and 0.6 GPa 


impact stresses, the rise times increase with propagation 


distance out to about 8 mm. It appears that the waves 


become steady after 7 to 8 mm of propagation distance at 


these impact stresses. At 1.2 GPa the rise time increases 


from 60 to about 100 ns between 2 and 4 mm. 


It is worth noting that the plastic wave in HMX single 


crystals 3 mm thick impacted at 1.4 GPa has a rise time 


of about 122 ns. In pentaerythritol tetranitrate (PETN) 


crystals unsteady plastic waves were observed.19 There- 


fore, some of the dispersed, unsteady wave behavior seen 


in PBX 9501 may be inherent in the HMX. 


The behavior of decreasing rise time with increasing im- 
pact stress has been noted before as the typical response 
to shocks. Swegle and Grady noted that for metals and 
oxides that the peak strain rate in the shock was consis- 
tent with a power-law relation between the peak strain 
rate and the shock stress;20 the peak strain rate for the 
total strain varies as the fourth power of the shock stress. 
They obtained the peak strain rate in the shock proﬁle by 
numerically diﬀerentiating the measured proﬁle. 


For comparing the PBX 9501 data to this relation, the 
strain rate in the wave front was taken as the total strain 
divided by the 10 to 90% rise time at sample thicknesses 
where the wave had become steady. 
The comparison is 
shown in Fig. 25. The line is 0.0395 times the strain rate 
to the one-fourth power. The data for PBX 9501 agree 
with the power law except at 0.3 GPa. Swegle and Grady 
had also noted this sort of deviation when the shock stress 


is not much larger than the elastic limit. 


I. 
Wave Speed as a Function of Particle Speed 


through the Unsteady Wave 


A steady wave is one whose proﬁle does not change as it 


propagates between two Lagrangian positions in the mate- 
rial. From the plot for rise time vs wave propagation dis- 


tance (Fig. 24), the increase in rise time with propagation 


distance shows that the wave is not steady at low impact 
stresses and short distances of propagation. The varying 


proﬁle can be expressed as the variation in wave velocity 
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Table 5: Rise Times for PBX 9501 Waves 


Shot 
Thick- 
Rise 
num. 
ness 
Px 
timea 


(mm) 
(GPa) 
(ns) 


1116 
1.322 
0.303 
400 


” 
2.318 
0.303 
– 


” 
3.298 
0.303 
740 


” 
4.27 
0.303 
900 


” 
5.274 
0.303 
1080 
1070 
5.415 
0.345 
1030 
1071 
5.405 
0.307 
1080 
1072 
7.993 
0.310 
1330 
1061 
7.993 
0.341 
1330 
1049 
1.394 
0.560 
228 
” 
2.412 
0.560 
266 
” 
3.405 
0.560 
372 
” 
4.374 
0.560 
406 
” 
5.370 
0.560 
496 
1117 
5.39 
0.582 
502 
1073 
7.976 
0.579 
660 
1060 
10.018 
0.598 
658 
1058 
1.37 
1.18 
62 
” 
2.36 
1.18 
72 
” 
3.36 
1.18 
97 
” 
4.33 
1.18 
110 


” 
5.33 
1.18 
106 
1054 
5.04 
1.20 
86 


1053 
5.016 
1.56 
28 


1062 
4.982 
1.61 
18 
1124 
10.015 
1.57 
50 


1125 
4.956 
1.80 
35 


1081 
4.538 
1.96 
17 


aThe rise times were measured from 10 to 90% of the total 


particle velocity increase. 
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Figure 24: Compressive wave rise time vs wave 
propagation distance for several impact stress lev- 
els in PBX 9501. 
The black exes are near 1.6 
GPa, the orange circle near 1.8 GPa, and the 
green circle near 2 GPa. 
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Figure 25: Comparison of the total strain rates 
of PBX 9501 for steady waves vs impact stress to 
the fourth power law. 
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with particle velocity through the proﬁle. The higher par- 


ticle velocity components are traveling more slowly than 


the low velocity components; so the wave spreads. Two 


examples for experiments at 303 and 560 MPa using data 


from wave proﬁles in Figs. 6 and 8 are given in Fig. 26. 


The particle velocity proﬁles at given Lagrangian gauge 


positions can be used to ﬁnd the phase velocities21 


Cu = ∂h 


∂t 


u 
. 
(7) 


For compactness of notation and consistency with previous 


work, u is used to denote longitudinal particle velocity and 


σ to denote longitudinal stress in the direction of wave 
propagation. Using the phase velocity Cu, the volume and 
strain through the wave, along the loading path, can be 
obtained using 


dV = − 
du 


ρ0Cu(u) . 
(8) 


Independent stress-time (σ−t) proﬁles are needed in order 
to obtain stress-strain curves at a given Lagrangian posi- 
tion h. The phase velocities at constant stress are deﬁned 
as 


Cσ = ∂h 


∂t 


σ 
. 
(9) 


The phase velocities Cσ and Cu are equal only for simple 
waves, steady waves or discontinuities. 
If they are ap- 
proximately equal, then one may substitute Cu for Cσ in 
Eqn. 10 


dσ = ρ0Cσdu 
(10) 


to calculate stresses through the wave. 


In order to obtain the phase velocities Cu, the arrival 


times at each of the ﬁve gauges one mm apart for a given 
particle velocity were recorded. In general the ﬁve data 


points appeared to follow a straight line within the preci- 


sion of the data. A linear ﬁt to the data yielded a wave 


velocity at that particle velocity. A constant value of Cu 
through the ﬁve Lagrangian positions (straight character- 
istics) deﬁnes a simple wave, in this case a nonisentropic 


simple wave.22 Therefore, the integrations for calculating 


V and σ with time at a given Lagrangian position can be 


carried out. 


At 303 MPa the separation into elastic and inelastic wave 


regions is evident in Fig. 26. The elastic range below 0.025 


mm/µs yields a ﬁt of Cu(u) = 2.81 − 5.94u. There is some 


wave dispersion even for the elastic wave. The inelastic 


range gives a ﬁt of Cu(u) = 3.16−19.5u. A similar analysis 


can be applied to the release waves. 


At 560 MPa the resultant wave velocities as a function of 


the particle velocity through the wave appear to fall on a 


straight line. A ﬁt to the data gives Cu(u) = 2.92 − 4.96u. 


Here the dispersion coeﬃcient 4.96 is close to the elastic 


coeﬃcient 5.94 at 303 MPa. The Cu(u) data at 560 MPa 


agree well with data obtained by Olinger and Hopson using 


manganin gauges.11 


While there is no strong nonlinearity in the wave arrival 
times at a given particle velocity between 1.3 and 5.3 mm, 
the waves appear to become steady by 8 mm of wave prop- 
agation according to Fig. 24. For a steady wave Cu(u) is a 
constant for all u; all parts of the wave are moving at the 
same velocity. In some way the transition from unsteady 
to steady must be made by all parts of the wave coming to 
a common velocity. In order to do a precise test for steadi- 
ness, it is necessary to perform an experiment to measure 
the in-material proﬁle at two Lagrangian positions in the 
steady region. If the wave is steady the transit time be- 
tween the two positions will be identical for all particle 
velocities. 
All states in the wave will lie on a Rayleigh 
line. One way for steadiness to be achieved is for the lower 
parts of the wave to slow down and the upper parts of the 
wave to speed up. 


J. 
Discussion of Experimental Errors 


For the 303 MPa experiment just discussed, the inter- 


cept at zero particle velocity is 2.81 mm/µs. This is be- 
low the typical value for the initial elastic disturbance of 


about 2.90 mm/µs by about 3%. It is low probably due 


to measurement error. It became apparent that it is sig- 


niﬁcant that each successive gauge that is 1 mm deeper in 


the sample is also staggered sideways by 1.7 mm (Fig. 1). 
This makes the experiment like a wedge test in that one is 


really measuring the phase velocity along the oblique face. 


This makes the measurement of wave velocity susceptible 
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Figure 26: Compressive wave phase velocity Cu 
as a function of particle velocity through the un- 
steady wave for two impact stresses. 


to wave tilt errors, especially at low impactor velocities as 
were used to get to low impact stresses. 


The tilt of the wave magniﬁes the tilt angle between the 
impactor and the sample surface by the ratio of the wave 
velocity to the impactor velocity. If the impactor velocity 
is 0.13 mm/µs and the wave speed is 2.6 mm/µs as is the 
case for 300 MPa symmetric impact experiments, the ratio 
is 20. The impactor tilt can be 1 mrad, where 1 deg=17 
mrad, although we attempt to keep it to a fraction of a 
milliradian. For 1 mrad of impact tilt the wave tilt is 20 
mrad or 1.2 deg. The ratio of the actual velocity U to the 
measured velocity UM is given by 


U 
UM 
= 
1 − sinα 


sinθ 


. 
(11) 


Here α is the impact tilt angle and θ is the angle that the 


gauge package is emplaced relative to the impact plane. 


For α of 20 mrad and a θ of 30◦ the ratio is 0.959. This 


could explain the error in wave velocity in experiment 1116 


at 300 MPa. It may also explain the high velocities ob- 


tained in experiment 1042. 


A number of experiments were performed with two in- 


plane stirrup gauges, one on the impact face and one em- 


bedded in the PBX 9501 at some depth. These measure 


impact time at the sample center directly and are not sus- 


ceptible to phase velocity errors due to tilt, as are the 


oblique embedded magnetic gauges, or the problems with 


PZT (lead zirconate titanate) pins at low velocities. There- 


fore, they should have the best precision for measuring 


transit times and wave speeds. 


In VISAR experiments the transit time was measured 


by the time diﬀerence between impact time and onset of 


the VISAR signal at the explosive/PMMA interface. The 


time of impact was measured by two pins on opposite sides 


of a diameter 54 mm apart. The average of the times mea- 


sured by the two pins should give the impact time at the 


sample center. In the experiments using stirrup gauges, 


1070 to 1073, the wave velocity was computed using the 
time diﬀerence between the impact gauge and the embed- 
ded gauge. 
In addition two PZT pins were used at the 
impact face in order to check the relative accuracy of the 
two methods. The samples were 5.4 to 8 mm thick. Im- 
pact velocities for the four experiments were 0.133, 0.134, 
0.149. and 0.148 mm/µs. The method using two stirrup 
gauges should be quite accurate. The PZT pin times have 
to be corrected for any deviation in position from the im- 
pact plane of the sample. 
The time uncertainty due to 
mechanical imprecision is magniﬁed at low impactor ve- 
locities. The average diﬀerence in transit times between 
the two methods of measurement was 0.11 µs. This gives 
diﬀerences in wave velocity of about 5% and 0.14 mm/µs 
for a wave speed of 2.6 mm/µs in the thinnest sample. 
We have started using four pins so that some statistical 
averaging can be performed to improve the precision. 


K. 
Spall Experiments on PBX 9501 at 600 MPa 


An attempt was made to measure the spall strength of 


PBX 9501. This is the strength in tension at high strain 
rates. A tensile state is achieved by the collision of rar- 


efactions in the sample after the initial compressive wave 


has passed through. For these experiments the sample rear 


surface is a free surface so that the compressive wave is re- 


ﬂected from it as a rarefaction to zero longitudinal stress. 
This rarefaction combines with the rarefaction from the 


back of the z-cut quartz impactor to generate tension. 


Experiments 1082 and 1084 were performed in this way. 
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The mirrors on the rear surface of the sample were 12 µm 


thick aluminum foils. The resultant proﬁle for experiment 


1084 is shown in Fig. 27. A typical spall signal shows a 


reduction in particle velocity due to the rarefaction origi- 


nating at the free surface. When spall occurs in the inte- 


rior the tensile stresses drop to zero there. When this is 


communicated to the free surface, an increase in particle 


velocity is observed. The magnitude of this rebound is a 


measure of the spall strength.23 There is no detectable re- 


bound in evidence in the experiments in PBX 9501. The 


inference is that the spall strength is zero within resolu- 


tion. Apparently the 600 MPa compressive wave commin- 


utes the explosive suﬃciently that it no longer has any 


strength; it acts like sand. 


At low strain rates in Brazilian tests the tensile strength 
has been measured at 0.78 ± 0.1 MPa24 at a failure strain 
of 0.0039. In Ref. 5 a value of 2.2 MPa is given for ultimate 
tensile strength. This tells us that the explosive is generally 
weak in tension. These values are much lower than the 
compressive strength of PBX 9501. 


Values this small might not have been resolvable in the 
spall experiments. One MPa is less than 0.2% of the com- 
pressive stress in the spall experiments. One can conclude 
that there is no large increase in tensile strength due to 
increased strain rate as there is in compressive tests. 


Spall behavior of ﬁlled elastomers has been studied by 
Kanel’ and coworkers.25 The samples had 75 wt % ﬁller. 
A small but distinct spall signal was noted. The particle 
velocity continued to decrease with time after the spall 
signal. This was interpreted to mean that the elastomer 


maintained continuity and continued to stretch after the 


initial void nucleation. Coarse ﬁller (160-200 µm particles) 


caused slower rise times and less deceleration than ﬁne ﬁller 


(20-50 µm particles). 


L. 
Additional Wave Proﬁles in PBX 9501 


For completeness, proﬁles from two embedded magnetic 


gauge, symmetric impact experiments are presented here. 


The ﬁrst experiment of this type near 600 MPa was per- 
formed on PBX 9501 that diﬀered from the explosive used 


in the remaining experiments. 
It was from another lot, 
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Figure 27: VISAR free-surface record of a spall 
experiment in PBX 9501 at 0.55 GPa, experiment 
1084. The sample was 3.01 mm thick. The z-cut 
quartz impactor velocity was 0.153 mm/µs. The 
impactor was 3.028 mm thick. 


was ram-pressed rather than hydrostatically pressed, and 
had a density of 1.835 g/cm3 as opposed to 1.827 ± 0.001 
g/cm3. It has shorter rise times and a higher elastic limit 
(Fig. 28). 


The ﬁrst experiment at 300 MPa was experiment 1059. 
The ﬁnal states measured were about 6% low (Fig. 29). 
The marked two-wave structure that it displayed was not 
repeated in subsequent experiments at this impact stress 
level. 
However, it made clear that an elastic wave and 
elastic limit existed. This was more subtle in subsequent 
experiments. 
One reason for the low precision in mag- 


nitude and proﬁle in this experiment is that the mating 


surfaces of the explosive where the gauge is emplaced were 
poorly machined. The surfaces were rough and not as ﬂat 


as is desirable. 


VII. 
EXPERIMENTS ON THE 


COMPONENTS OF PBX 9501 


Numerical computation has advanced to a stage that 


wave propagation in a composite can be computed on a 


heterogeneous matrix of the components. Therefore, it is 


worthwhile to measure the plate impact response of the 
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Figure 28: In-material MMV gauge records in 
PBX 9501 at 598 MPa in a symmetric impact ex- 
periment, 1042. The PBX 9501 was ram-pressed 
at TA-9 from Holston lot 85G730-007. Gauge po- 
sitions are 1.38 to 5.32 mm in from impact face. 
Half the impact velocity is 116 m/s. 
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Figure 29: In-material MMV gauge records in 
PBX 9501 at about 300 MPa in a symmetric im- 
pact experiment, 1059. Particle velocities are low, 
and the proﬁles appear distorted. However, they 
show the two-wave structure clearly. Gauge po- 
sitions are 1.43 to 5.38 mm in from impact face. 
Half the impact velocity is 63.4 m/s. Sample den- 
sity was 1.828 g/cm3. 


individual components. PBX 9501 is 95 wt % HMX in a 


bimodal distribution and 5 wt % binder. The binder is 


half estane and half nitroplasticizers. A small amount of 


antioxidant is added also. 


A. 
Hugoniot Data and Wave Proﬁles for PBX 


9501 Binder and Estane 


Impact experiments were performed on the binder and 


estane from 200 to 1000 MPa. The binder is half estane, a 


polyester polyurethane elastomer, by weight. The results 


for Hugoniot states are summarized in Table 6. 
In the 


binder experiment 1085 at 1000 MPa an accurate wave 


transit time was not obtained, since only one of the two 
PZT impact pins was recorded. 
However, a good wave 
proﬁle was measured. 


Hugoniot data for estane from 1.1 to 18 GPa is given in 
Ref. 9. The ﬁt of shock velocity vs particle velocity for 
that data is US = 2.32 + 1.70up. The intercept is signif- 
icantly higher than the bulk sound speed of 1.75 mm/µs 
(Table 2) computed from our measurements of longitudi- 
nal and shear sound speeds. In a simple material the bulk 
sound speed should agree with the intercept of the bulk 
Hugoniot. The previous data and the new data are plot- 
ted in Fig. 30. The new data clariﬁes the situation. They 
show that the low pressure Hugoniot has a diﬀerent slope 
than the high pressure Hugoniot, and that the low pressure 
data agrees with the bulk sound speed as intercept. This 
clariﬁes a discrepancy noted in many polymers.26 


The estane data is shown as pressure vs volume in 


Fig. 31. It is quite compressible up to about 15% and then 


stiﬀens. Apparently the free volume is taken out so that 
interchain repulsions dominate at higher pressure. 
This 


may be aﬀected by the degree of crosslinking. The Hugo- 


niot data has been included in a thermodynamic Helmholtz 
potential by Ralph Menikoﬀ.27 The low pressure data was 


described by treating estane as having 1.5% porosity in a 


Px−α model. This ‘porosity’ may be thought of as existing 


at the interchain level. 


The Hugoniot data for the binder are displayed in 


Fig. 30. The points extrapolate back to the bulk sound 


speed in a fashion similar to estane. It is to be expected 
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Table 6: Results of VISAR Experiments on PBX 9501 Binder and Estane 


Shot 
Sample 
Impactor, 
Thick- 
num. 
typea 
d 
ness 
uI 
US 
upb 
Pxb 
ϵxb 
T 


(mm) 
(mm) 
(mm/µs) 
(mm/µs) 
(mm/µs) 
(GPa) 
(◦C) 


1088 
binder 
TPX, 5.33 
6.16 
0.196 
1.94 
0.080 
0.198 
0.041 
17.0 


1066 
binder 
z-qtz, 6.37 
4.49 
0.179 
2.13 
0.154 
0.420 
0.072 
20.2 


1085 
binder 
z-qtz, 6.37 
5.48 
0.373 
– 
– 
– 
– 
– 


1087 
estane 
TPX, 5.55 
10.018 
0.1885 
2.11 
0.076 
0.191 
0.036 
17.0 


1065 
estane 
z-qtz, 6.37 
5.07 
0.145 
2.31 
0.124 
0.345 
0.054 
19.0 


1080 
estane 
z-qtz, 6.37 
5.01 
0.347 
2.75 
0.291 
0.959 
0.106 
19.8 


aThe binder samples were made at S-Site by pouring thin layers of lacquer into muﬃn tins and letting the solvent 
evaporate between pours. The binder material was assayed at 51.8 wt % estane and 48.2 wt % nitroplasticizer. By 
GPLC (gel-permeation liquid chromatography) the peak and average molecular weight of the estane in the binder were 
119,000 and 153,000 g/mol, respectively. The estane samples were all hot-pressed under vacuum to about 1 ksi at S-Site 
from estane ﬂakes from B. F. Goodrich Company. The pressed billets were 6 inches in diameter and 2 inches thick. The 
estane lot was 6812100SK13. The average and peak molecular weights were about 73,000 g/mol. 


bValues were calculated from the intersection of the Rayleigh line and the appropriate reﬂected impactor Hugoniot using 
the wave velocity in the sample at half maximum particle velocity, the measured impactor velocity and the known 
impactor Hugoniot. 
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Figure 30: Shock velocity vs particle velocity 
Hugoniot data for estane and PBX 9501 binder. 
The blue points are the new data and measured 
bulk sound speed. The black triangles are from 
Marsh’s compendium. 
The red exes are bulk 
sound speed and Hugoniot points for the binder. 


0 


5 


10 


15 


20 


0.6 
0.7 
0.8 
0.9 
1.0 
1.1 


Stress (GPa) 


Relative Volume 


Figure 31: 
Pressure-volume Hugoniot for es- 
tane. 
The new data are the blue points. 
The 
others are from Marsh’s compendium. 
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that the US vs up Hugoniot will have a lower slope at 


higher pressures as is the case for most polymers including 


estane.26 


The wave proﬁles from the VISAR experiments on the 


binder are displayed in Fig. 32. 
The window material, 


PMMA, has a higher shock impedance than the binder or 


estane. Therefore, shocks are reﬂected from the interface 


back into the sample, instead of the rarefactions experi- 


enced by PBX 9501 samples. The impact stresses for the 


three experiments are 198, 420 and 1040 MPa. 


Using MACRAME and the measured Hugoniots the cal- 


culated interface particle velocities for the three experi- 


ments are 67.5, 135, and 293 m/s. The measured states be- 
hind the shock are 72, 123 and 291 m/s, respectively. The 
intermediate experiment has a sharper wave front than the 
other two. The other two have more rounding before the 
ﬁnal state. These two are probably more representative of 
the material. 


The 420 MPa experiment had a VISAR fringe constant 
of 176 m/s per fringe, and the 198 MPa experiment had 
a fringe constant of 96 m/s per fringe. Thus both proﬁles 
were obtained from about 3/4 of a fringe. The VISAR par- 
ticle velocity accuracy depends on the number of fringes. 
With less than a fringe the accuracy is not good. 
The 
highest stress proﬁle has three fringes in the shock rise and 
should be accurate to about 1%. The gradual ramping of 
the low stress experiment and the late bumps on the ex- 
periment at the intermediate stress are probably artifacts 


of the experimental imprecision with less than a fringe of 
data. 


There is 250 MHz noise on the fringes due to shot noise 


in the photomultiplier tubes. For some experiments the 
fringe data was smoothed 5 points at a time to reduce 


the high frequency noise on the signal. The analysis pro- 


gram written by Will Hemsing avoids smoothing the sharp 


jumps. 


Similarly, the estane proﬁles from VISAR experiments 
are shown in Fig. 33. The impact stresses determined from 


the measured shock speeds and the impactor velocity and 


Hugoniot are 191, 345, and 959 MPa. Using MACRAME 


and the measured Hugoniots the calculated interface par- 
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Figure 32: VISAR interface records for PBX 
9501 binder at the interface with PMMA at 198, 
420, and 1040 MPa impact stresses. For the ex- 
periment at 1040 MPa the wave arrival time is 
estimated. 


ticle velocities for the three experiments are 65, 111, and 
271 m/s. The measured states behind the shock are 66, 
105 and 272 m/s, respectively. 


As with the binder the two lower stress experiments had 
less than one fringe of data. The apparent oscillations are 
probably not due to the inherent impact response of es- 
tane. Besides the problem of experimental imprecision the 
samples were soft and diﬃcult to machine to good toler- 


ances of ﬂatness and thickness uniformity. In assembling 


the target, epoxy was used to ﬁll in any gaps due to devi- 


ations from ﬂatness in the sample. Generally the samples 


were thickest in the center so that epoxy thickness was 
minimized at the laser spot on the mirror. The resultant 


eﬀect of these sample imperfections on particle velocity 


measurement is magniﬁed at low impact velocities as used 
for the experiments at 191 and 345 MPa. 


The experiment at 959 MPa should be accurate. The 
ramping rise behind the shock is quite marked, but the 


ﬁnal state is in good agreement with a MACRAME calcu- 


lation using measured Hugoniots. 
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Figure 33: Particle velocity proﬁles for estane 
at the interface with PMMA. Impact stresses are 
190, 340, and 959 MPa. 


B. 
Experiments on HMX Single Crystals 


Two VISAR experiments with PMMA windows were 
performed in order to measure the Hugoniot elastic limit of 
β-HMX crystals. The 3 mm thick crystals were mounted 
on 5 mm thick Kel-F discs and surrounded with an 
impedance matching mixture of epoxy and silica. The im- 
pact stress was about 1.4 GPa. The two orientations stud- 
ied had faces parallel to the {110} and {011} crystal planes 
in P21/n space group. HMX β phase is monoclinic. Crys- 
tallographers use P21/c; crystal growers use P21/n, while 
the slip systems of anthracene and other molecular crystals 
of the same space group are given in P21/a. The revised in- 


ternational crystallography tables28 list space group P21/c 


with the other two as alternate choices of a coordinate sys- 


tem. 


The orientations were determined by measuring inter- 
facet angles with a protractor. 
By comparing a set of 


angles with computed ones, it was possible to determine 


the facets. Slabs were cut parallel to a facet with a di- 


amond saw to prepare a sample for impact experiments. 


The {110} crystal was 3.210 mm thick, and the {011} 
crystal was 2.997 mm thick. The proﬁles at the interface 


with a PMMA window from the VISAR experiments are 


shown in Fig. 34. The {110} orientation has an elastic- 


plastic proﬁle similar to that seen in simple materials 


with some work hardening between the elastic and plastic 


waves. 
The {011} orientation has a stress-relaxing pro- 


ﬁle similar to that observed in LiF, Armco iron, and other 


materials.29, 30, 31 The elastic shock strength decays until 


a steady value is reached. 


Howard Cady performed early studies on quasi-static 


deformation of HMX. Twinning is an important deforma- 


tion mechanism in β-HMX. The twin plane is (101).32 The 


twin direction has been postulated to be [10¯1].33 Other de- 


formation mechanisms may also contribute. The fracture 


plane is on planes of {011} type.32 Slip planes (100) and 


(101) have been established with possible Burgers vectors 
of [100] and [10¯1].33 The relationship of the impact plane 
to the deformation systems (plane and direction) will de- 
termine the elastic-plastic behavior. Here, the word ‘plas- 
tic’ is used loosely to cover the inelastic deformation. In 
an isotropic material the maximum resolved shear stress 
in uniaxial strain is at 45◦ to the impact plane. The defor- 
mation system that has plane and direction closest to this 
will likely take most of the inelastic deformation, once the 
critical shear stress is exceeded. If the deformation mecha- 
nism has signiﬁcant kinetics to it, stress relaxing behavior 
will be observed. 
Apparently diﬀerent mechanisms and 
kinetics dominate in the two orientations studied. 


Elastic precursor amplitudes were determined by using 
Rayleigh lines for the crystal elastic release isentrope to 
get from the interface state to the state in the crystal. 


For the {110} crystal the velocity of the ﬁrst disturbance 
was 4.25 mm/µs, and the velocity of the elastic wave at 


half maximum particle velocity was 4.16 mm/µs. This is 


substantially larger than the longitudinal sound speed (Ta- 
ble 2). This implies that the crystal is nonlinearly elastic; it 


stiﬀens in response to increasing impact stress. Similar be- 


havior was observed in PETN crystals.34 The elastic state 
in HMX was determined to be 460 MPa and 58 m/s. The 


plastic wave had a Lagrangian wave speed of 3.28 mm/µs 


and a rise time from 10 to 90% of 123 ns. 
The plastic 


waves in PETN were unsteady.19 


For the {011} crystal the speed of the initial disturbance 
was 4.07 mm/µs and 4.04 mm/µs for half maximum on 
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Figure 34: 
Particle velocity proﬁles at the 
HMX/PMMA interface from VISAR experiments 
for crystals with [110] and [011] faces. The elastic 
precursor wave is evident. 


the elastic wave. In comparison, the measured longitudi- 
nal sound speed is 3.72 mm/µs. The elastic state in this 
orientation was computed to be 490 MPa and 64 m/s. At 
the HMX/PMMA interface the elastic state relaxes from 
86 m/s down to 62 m/s in about 70 ns. The plastic wave 
has a rise time of 123 ns, the same value as the other ori- 
entation. The measured Lagrangian plastic wave was 3.21 
mm/µs speed, very close to the value for the {110} orienta- 
tion. The Lagrangian wave speed converts to an Eulerian 
wave speed of 3.22 mm/µs. So the two orientations have 
similar but not identical elastic-plastic responses to im- 


pact. More experimental work is needed to delineate the 
details of the behavior. 
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